Orbital-transverse density-wave instabilities in iron-based superconductors 



O 

(N 



(N 



I 

C 

o 
o 



> 

o 

o 
o 



Zi-Jian Yao, 1 Jian-Xin Li, 2 Q. Han, 3,1 and Z. D. Wang 1 

'Department of Physics and Center of Theoretical and Computational Physics, 
The University of Hong Kong, Pokfulam Road, Hong Kong, China 
* National Laboratory of Solid State Microstructures and Department of Physics, Nanjing University, Nanjing 210093, China 
s Department of Physics, Renmin University of China, Beijing, China 
(Dated: July 27, 2010) 

Besides the conventional spin-density-wave (SDW) state, a new kind of orbital-transverse density- 
wave (OTDW) state is shown to exist generally in multi-orbital systems. We demonstrate that the 
orbital character of Fermi surface nesting plays an important role in density responses. The relation- 
ship between antiferromagnetism and structural phase transition in LaFeAsO (1111) and BaFe2As2 
(122) compounds of iron-based superconductors may be understood in terms of the interplay be- 
tween the SDW and OTDW with a five-orbital Hamiltonian. We propose that the essential difference 
between 1111 and 122 compounds is crucially determined by the presence of the two-dimensional 
d^y-like Fermi surface around (0,0) being only in 1111 parent compounds. 

PACS numbers: 75.25.Dk, 75.30.Fv, 74.70. Kn 



I. INTRODUCTION 

Over the last two years, research on iron-based super- 
conductors has been an exciting topic that attracts in- 
tensely experimental [3-[E] and theoretical [^-[Tl| inves- 
tigations. Although it has been well established exper- 
imentally that this family of compounds exhibit several 
phase transitions including the structural phase transi- 
tion 
tion 



.12 4 141 , the antiferromagnetic (AF) phase transi- 
3 15j |. and the superconducting phase transition, 



the mechanism of these phase transitions remains highly 
controversial. However, there has been a consensus on 
the basic Fermi surface (FS) topology - hole pockets cen- 
tered at (0, 0) and electron pockets centered at (tt, tt) flB - 
[l8j . From a viewpoint of itinerant antiferromagnetism, 
the hole and electron pockets are assumed to be nested 
nearly perfectly in the parent compound. The spin- 
density wave (SDW) state is stabilized due to the ex- 
istence of the on-site Coulomb repulsions. Upon doping, 
the long range AF order is destroyed and short range AF 
spin fluctuations are developed, which is responsible for 
the high-temperature superconductivity (l9l [20j . 

The multi-orbital nature of iron-based superconduc- 
tors is believed to play a prominent role in the supercon- 
ductivity, which signifies the importance of extracting the 
distinct physics that emerges from multi-orbital effects in 
this new family of materials. In this paper, we reveal that 
a new kind of unconventional orbital-transverse density- 
wave state exists generally in multi-orbital systems with 
certain orbital configuration of FS nesting. This type 
of density wave stems from the rotation asymmetry of 
Hamiltonian in the orbital space, which is reflected in the 
multi-orbital FS nesting geometry, and has an intriguing 
impact on the electron charge and spin density responses. 
In connection to the iron-based superconductors, our cal- 
culations show that the orbital-transverse density wave 
is a competing order with conventional SDW. 

This paper is organized as follows. In Section [Til 
we first consider three prototypes of FS nesting in the 



multi-orbital systems. Then we define the density oper- 
ators with two orbital degrees of freedom, and show how 
their responses are affected by the orbital configuration 
of FS nesting. In Section IIII1 we show that the orbital- 
transverse density wave is nearly degenerated with the 
SDW by using a realistic five-orbital Hamiltonian of the 
iron-based superconductors. The interplay between the 
structure distortion and antiferromagnetism is discussed 
in the context of competing density- wave ground states. 
Finally, several remarks are drawn as the summary in 
Section HVl 



II. DENSITY WAVES WITH ORBITAL 
DEGREES OF FREEDOM 

A. Multi-orbital Hamiltonian 

A multi-orbital Hamiltonian incorporating the on-site 
intra- and inter-orbital Coulomb interaction reads 



-iota ,tj %a.<y 



(i) 



where H = S ic ^ CT ii Q j/3q aff c^ CT , Ua,jp is the hopping 
term between orbital a of site i and orbital /3 of site j, and 
rii aa the electron number operator, U and U' is the on- 
site intra- and inter-orbital Coulomb repulsion, respec- 
tively. For simplicity, U is assumed to be equal to U' in 
the present study and the terms of Hund's rule coupling 
are ignored. In order to illustrate various types of den- 
sity responses of the multi-orbital system more clearly, 
we start from a simple two-orbital model in the square 
lattice with a hole pocket lying at (0,0) and an electron 
pocket lying at (tt, 0). The dispersion and chemical po- 
tential ensure that the electron pocket and hole pocket 
are circular and nested perfectly. This fermiology is com- 
patible with many possible orbital configurations, which 
have different orbital-weights of the FS. In Fig. 1 we 




(a) 



(d) 



FIG. 1: (Color online) Schematic plot of the Fermi surfaces 
with different orbital configurations. Blue and red indicate 
dominant weights of orbital I and 2, respectively. 



show three prototypes of them. As will see below, the 
density response behaviors of them are quite different. 



B. Definition of density quantities 

The introduction of orbital degrees of freedom enables 
us to define extra physical quantities of density whose 
translational symmetry could be broken when the sys- 
tem is phase transferred into the density-wave ground 
state. Simiar to the spin-i degrees of freedom, we repre- 
sent the two-orbital degrees of freedom by an i pseu- 
dospin. The density operators in such a spin-orbital 
space can be written as D(rj) = (j>\T<j)i with T = r ® a 
and <f>\ — [cj^, cj^, c| 2 ^, c^ 2 jj the 4-component spinor 
(t° = a = I and r 1-3 = er 1-3 the Pauli matrices, 
where t is defined in the orbital space and a is de- 
fined in the spin space). To study the instability in the 
particle-hole channels, the sixteen particle-hole operators 
<jr (t 0-3 (£)<T a ~ 3 )(f> can be reduced into six different density 
channels (three of spin and three of charge) due to the 
SU(2) symmetry of Hamiltonian {1} in the spin space, 
which are listed in Table H 

The six channels listed in Table Q] could be further 
grouped into four types of density waves, if we are 
not interesed in the direction of the density orbital- 
polarization, namely the total spin-density wave, the 
total charge-density wave, the orbital-polarized spin- 
density wave, and the orbital-polarized charge-density 
wave. Basically the orbital-polarized density operator 
represents the difference of density between two orthog- 
onal orbitals. The longitudinal or transverse polariza- 
tion of orbital determines which two orthogonal orbitals 
are density-ordered, for instance, the longitudinal polar- 
ization corresponds to a density-ordered state between 
orbital 1 and 2 while the transverse polarization corre- 
sponds to a density-ordered state between two orthogonal 
mixtures of orbitals I and 2. The distinction between 
orbital-longitudinal and orbital-transverse parts of the 
susceptibilities is consequent to the rotation asyemmetry 
of the Hamltonian (fT]) in the orbital space. In Fig. 2 
we plot schematically the ground states of these density- 
wave states in real space. 



AA 


tt 


AA 


tt 










•J 


W-..J 




•J 










II 


II 


II 


II 


1 1 




1 1 








(b) 






(e) 






At 


tA 


At 


tA 










• 


•o 

II 


•o 
II 


•o 

II 


mo 

r 


•3 

'1 


•o 

I 1 


mo 
I | 






(c) 






(i) 






At 


tA 


At 


tA 










II 

i-1 


II 

i 


II 

1+1 


II 

i+2 


r 

i-i 


'1 

i 


I 1 

i+1 


I | 

i+2 



FIG. 2: Schematic plot of density waves with wave vector 
Q = (tt, 0), (a)-(f) correspond to the density-wave states 
with density operators Do, 3, ^3,3, 3> Do t o, ^3,0, and 

-D_(+) i o. Black (blank) circle represents orbital 1 (2). Dark 
grey and grey circles represent two orthogonal mixtures of or- 
bitals 1 and 2, e.g., -L(|l) + |2)) for dark grey and ^(]1)-|2)) 
for grey. The arrow indicates spin density and the vertical line 
indicates the charge density, i, i + 1, ... denote the lattice sites 
and each site represents a reduced unit cell. Only modulated 
direction is shown. 



C. Instability criterions 

The particle-hole instabilities in the presence of 
Coulomb interactions may be examined within the 
random-phase approximation (RPA). The interaction 
matrices for the magnetic (U s ) and charge (U c ) chan- 
nels reads 




(2) 



where U sl = diag(J7, U), U£ n = U for m = n and 2U' 
otherwise, II s2 = diag(Z7', U') = ~U c2 . Defining the 
orbital indices of the bare susceptibility as XaSuvil) = 
-|E t ^(H?)^W, we get the RPA spin and 
charge susceptibility matrices 

x c (q)^[i + x°(q)uT 1 x°(q)- (3) 

The response functions listed in Table Q] are given by 
1 rfi 

XA, B {q,iu) = - J dTe WT (T T D A , B (q,T)D AtB (-q,0)) 
1 rf> 

X- +lB (qM = g / dTe™ r {T T D_ iB (q,T)D +tB (-q,Q)) 

(4) 

where A and B are numbers and D a,b = 4^ { t a®&b) ( I ) - A 
direct calculation gives the response functions of density 
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TABLE I: Density waves with two orbital degrees of freedom 



Ordering type 
of densities 


Orientation in 
orbital space 


Density operator 


xvesponse function 


Divergent 
blocks 


total spin 


n/a 


-Do, 3(^) 


Pilt — Pill + P»2t — Pi2l 


Kxfni + X2222 + X1122 + X2211) 


x sl 


total charge 


n/a 


Do,o(ri) 


Pilf + PiH + Pi2f + Pi2l 


|(xini + X2222 + X1122 + X2211) 


x cl 


orbital- 


longitudinal 




D 3 , 3 (ri) 


Put - Pm - (pin ~ Pin) 


2 (Xnn + X2222 — X1122 — X2211) 


x sl 


polarized spin 


transverse 


D_ 




c i2t c *it - cLi^n ( H - c -) 


X2121 


x s2 


orbital- 


longitudinal 




o 3 ,o(n) 


Put + Pm - (p<2t + ^224.) 


|(Xini + X2222 — X1122 — X2211) 


x cl 


polarized charge 


transverse 


D. 


(+),o0i) 


4 2 t c »it + C i24. C iU (H.c.) 


X2121 


r 2 



operators defined in Table Q] 



X0,3 


— 2~(x!lll " 


" X2222 + X1122 


+ X2211 


Xo,o 


= i^oCwii ' 


l~ X2222 + X1122 


+ X2211 


X3,3 


— 2^1111 " 


" X2222 — X1122 


— X2211 


X-+.3 


= X2121) 






X3,0 


= l^oCwu - 


" X2222 — X1122 


— X2211 


X-+,o 


= X2121- 







(5) 



To a good approximation, especially at Q — (tt, 0), the 
bare susceptibility matrix in the two-orbital case can be 
regarded as a block matrix x° = diag(x\ X 2 ), where 



■ _ , , V 2 = l Xl212 Xl221 

x I I ' X I 

\ X2211 X2222 / V X2112 X2121 





X1111 X1122 



As a result, the RPA spin and charge susceptibility ma- 
trices are also blocked, 



>S (c) 



£al(cl) 




(6) 



where 



x sl 



Ml v s2 - Mi v cl -Mi^_Ml m 
di' x d 2 ' x d 3 ,x d 4 - {,) 



Here M n are some 2x2 matrices and 

di = det[x sl ] = (1 - f7xim)(l - UX2222) - ^ 2 Xii22X°2ii, 
d 2 = det[ X s2 ] = (f-(/' Xl212 )(l-t/'x ( 2 , 121 )"C/' 2 Xi22iX2ii2, 

di = det[x cl ] = 1 + ^(xi'in + X2222) + 2(7' (X1122 + 2x°2ii) 

. ittI ATT'2\r 00 \ 

+(U -4(7 )(XimX2222 - X1122X2211) 

d A = d 2 (8) 

are the denominators that come from the inversion of 
matrices I - x a (l)U s and I + x°(q)U c in Eq. ©. Devi- 
ated from the single band RPA, there are two instability 
criterions for the spin channel, which are indicated by 
d\ — > and c?2 — > 0, respectively. When approaching the 
critical temperature, either the upper block x sl or the 



lower block x s2 may diverge. Comparing the response 
functions listed in Table Q] and Eq. 0, we can see that 
the divergence of the upper block or the lower block cor- 
responds to a transition to one or the other totally dif- 
ferent magnetic phase, i.e., the former indicates an SDW 
transition while the latter indicates an orbital-transverse 
polarized spin-density- wave transition. The same conclu- 
sion holds for the charge susceptibility matrix. We will 
refer the orbital-transverse density-wave ground state as 
OTDW in this paper hereafter. 

Now we turn to the real-space distribution of these den- 
sity waves. As discussed in Section IIIB| it is clear that 
the orbital-longitudinal operators represent the density 
differences between orbitals 1 and 2. On the other hand, 
the real-space distribution of the orbital-transverse den- 
sity wave can be revealed by performing a rotation in the 
orbital space: 



~ [ C i+t' C i+-|,' C i-t' C i-4.]' 



where 



C l±a ~ _^2 ( C ilo' = ' =C i2cr)- 

The charge and spin density differences between the two 
orbitals |+) and |— ) can be represented by 

D 3,3( r t) = ip\ (T3 ® a 3 )ip = S i+iZ - Si- <z 

= ct lt C i2 t - 4ll C iU + HC ' 
^3,0 ( r ») = i>i(T3®(ro)lp = Pi+ ~ Pi- 

= c| lt ci2t + 4x^21 + H.C, 

where Sf± z and pi± represent the spin and charge den- 
sity of orbital |±) at i site, respectively. Similarly, the 
response functions for density operators D' A B are 



X3 : 3 — 2 (Xl212 + X2121 + Xl221 + X2II2)' 

X3,0 = 2^1212 +X2121 +X1221 +X2112)- (9) 

Therefore, in OTDW state, the translational symmetry 
of the local density difference between the two orthogonal 
mixtures of orbitals 1 and 2 is broken. 



D. Nesting-induced OTDW: a weak-coupling 
analysis 

Roughly speaking, the density response of a Hubbard- 
like Hamiltonian is mainly determined by two factors: 
one is the band structure, and the other is the Coulomb 
interactions. We now look into the influence of different 
band structures shown in Fig. 1 on the density responses. 
The bare susceptibility of a general multi-orbital system 
is give by 



X 



VM) = Jf E <{k + q)ai{k)al*{k)a^{k + q) 

k,mn 

f[e n (k + q)]-f[e m (k)} 



e m (k) - e n (k + q) + ir)' 



(10) 



where a, /?, fi and v are the orbital indices, m and n are 
the band indices, and a is the orbital weight. For the 
case indicated in Fig. 1 (al), the matrix elements of bare 
susceptibility are simplified to 



f[e v {k + q)}- f[e,{k)\ 
V ^ ^{k)-e„(k + q)+i7 1 ' 



(11) 



The susceptibilities x'uuvu ^ v) vanish since there 
is no hybridization between the two orbitals. Straight- 
forwardly, Xi2i2( ( 3) an( l X2i2i(Q) are divergent due to 
FS nesting, while Xim an( l X2222 keep finite as they 
take no advantage from the nesting between electron and 
hole pockets. Bear in mind that although the cases of 
Fig. 1 (al) and (a2) have different orbital configura- 
tions, their FS nestings are both from the inter-orbital 
contribution [2l[ . Hence the above conclusion holds also 
for the band structure shown in (a2). For the case Fig. 
1 (b), from Eq. (ITU)) the divergent components of bare 
susceptibilities are X1111 an( l X§222- 

The above considerations focus on the band structure, 
but do not take into account the electron-electron interac- 
tions explicitly. Given the interaction vertices in Eq. ([2]) 
and the instability criterions in Eq. ([8]) , three spin-density 
waves (-Do, 3> -D33, and -D_( + ) 3 ), and one charge-density 
wave (£>_ (^0) could be established in Hamiltonian (fTJ. 
The I?o,o an d -D30 charge-density waves are prohibited 
due to repulsive on-site Coulomb interactions (note that 
c?3 is always greater than zero when ^Xim/2222 < 1 anc ^ 
the hybridization is neglect able ) . D0.3 and -D33 density 
waves occur simultaneously since they both come from 
the divergent of the upper block of X s (Q), which indicates 
the real-space distribution of spin density to be compat- 
ible with Fig. 2 (a) and (b). A typical pattern of spin 
density that is either -Do. 3 or D33 density wave is shown 
in Fig. 3, where the staggered spin density is orbital- 
polarized and dominantly appears in one orbital. On the 
other hand, the phase transition in -D_( + \ 0(3) channel 
comes from the divergent of the lower block of x s ^°HQ) 
which is coupled to the inter-orbital repulsion, with order 
parameter (cj^fc + Q)c l/ff ( (T )(fc)). In the present model, 



!♦ K !♦ K 



i-1 



i 



i+1 



i+2 



FIG. 3: Intra-orbital nesting induced spin-density-wave state 
that is realized in the iron-based superconductors, where the 
translational symmetries of {-Do, 3) and {1)3,3} are broken si- 
multaneously. All the symbols are the same as that in Fig. 
2. 



instabilities in D—i+\ 3 and channels are degener- 

ated due to U = U', while the actually established state 
depends on the additional interaction terms that break 
the degeneracy, e.g., the former is favored by Hund's rule 
coupling while the latter is favored by the presence of 
electron-phonon interactions. 

We now conclude that for configurations shown in Fig. 
1 (al) and (a2), the divergent susceptibility matrix is 
X s2 (Q) or X c2 (Q), indicating that the system undergos 
an OTDW transition. There is no long range order of the 
total local spin density, as is shown in Fig. 2 (c) and (f). 
For configuration Fig. 1 (b), the divergent susceptibility 
matrix is X s1 (Q), corresponding to a spin density- wave 
transition whose real-space distribution is shown in Fig. 
3. 



III. OTDW AND THE IRON-BASED 
SUPERCONDUCTORS 

The multi-orbital band structure of iron-based super- 
conductors provides an intermediate example between 
prototypes of Fig. 1 (a) and (b). To be more concrete, 
we here employ a realistic five-orbital Hamiltonian with 
on-site intra- and inter-orbital Coulomb repulsions. The 
involved orbitals are d xz , d yz , d x 2_ y 2, d xy , and d 3z 2_ r 2 
(defined in the reduced unit cell) which are labeled with 
1, 2, 3, 4, and 5 respectively. Then the orbital indices 
in Hamiltonian (fTJ) run from 1 to 5. We adopt the hop- 
ping parameters given by Graser et al [22| fitted from 
first principle calculated band structure by Cao et al [23| . 
The density- wave instabilities defined in Table Q] can be 
realized in the space spanned by the two most relevant 
orbitals. In our calculations, d xz /d yz (orbital 1/2) and 
d X y (orbital 4) are the two orbitals with largest intra- and 
inter-orbital susceptibilities. The instabilities of OTDW 
and SDW are found to be nearly degenerated, conse- 
quently the competition between them is sensitive to 
the detailed band structure. We consider this to be a 
key point to understand the essential difference between 
i?FeAsO (1111) and AFe 2 As 2 (122) compounds (i?=rare 
earth and ^4=Sr, Ca, Ba and K), which will be discussed 
in the following paragraphs. 
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The AF and structural transitions in the iron-based 
superconductors exhibit different properties in different 
systems: the static AF order develops after the struc- 
ture distortion in 1111 compounds [13], while in the 122 
compounds, the two transitions occur at the same tem- 
perature [lj, [H[ • The structure distortion has been 
theoretically proposed to be driven magnetically pBI - 
[27| or electrically (2|| [29| . Here we present an alterna- 
tive view. Note that all the density waves we discuss 
are nesting-driven and occur at the same wave-vector 
Q = (it, 0). The order parameters of these density waves 
and the tetragonal-to-orthorhombic distortion break the 
same symmetry. Therefore, we propose that either spin- 
polarized SDW (shown in Fig. 3) or OTDW could induce 
a structure distortion. When the transition temperature 
of SDW T dl is higher than that of OTDW T d2 , the AF 
transition is orbital-polarized as shown in Fig. 3 and is 
accompanied with a structure distortion which is mag- 
netically driven, and the critical temperature of struc- 
ture distortion T$ and AF transition Taf equal to Td\- 
The temperature gap between Tg and Taf appears when 
Td2 > T d i, which indicates T s = T d2 > T A f = T dl . We 
make two remarks here. First, D_( + ) i o and D_( + ) j 3 den- 
sity waves are two candidates for OTDW-induced lat- 
tice distortion. As an orbital-ordered state of charge, 
-D_( + ),o density wave could induce a structure distortion 
electrically by Jahn- Teller effect. The coupling between 
-D_( + ),o density wave and lattice distortion is, if exists, 
magnetic relevant and more subtle. It is out of scope 
of the present work to determine which of the two den- 
sity waves is actually established in the real materials 
of iron-based superconductors. Since the two states are 
degenerated in our model, we just generally related the 
structure distortion with OTDWs. Second, in this sce- 
nario, from high temperature to low temperature, the 
system undergos two phase transitions, both couple to a 
structure distortion. However, only one structure distor- 
tion is identified in experimental measurements. We note 
the lattice distortion is weak and driven by breaking of C4 
lattice rotation symmetry. Thus if the system is already 
in the orthorhombic structure phase, another density- 
wave phase transition with wave- vector (tt, 0) would not 
break an extra lattice symmetry. 

It is known from band calculations of 1111 compounds 
that along T-Z direction lies a d^-like two-dimensional 
(2D) FS (denoted as e-FS hereafter) [H- While for 122 
compounds, a 2D FS with d xy character is found to be 
absent when the pnictogen height is relaxed using total- 
energy minimization [3l|. In our model calculation we 
find the existence of e-FS enhances the inter-orbital sus- 
ceptibility X1414 substantially, while without e-FS, the 
susceptibility X1111 is i n general, at least sightly, stronger 
than that of X?4i4- The enhancement of X1414 by the 
appearance of a d xy -\ike FS is expected because of the 
strengthened inter-orbital nesting between e-FS and the 
4r.z/^2/z-dominant FS sections around M. In Fig. 4 we 
show dominant components (1111 and 1414) of spin sus- 
ceptibility matrix with (upper) and without (lower) the 
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FIG. 4: Dominant components of spin susceptibility matrix 
with (upper) and without (lower) the appearance of e-FS. The 
hopping parameters are adopted from |22(, U = U' = 1.38 and 
T = 0.19 for the upper panel and U = U' = 1.53 and T = 0.17 
for the lower panel. To produce the e-FS, the on-site energy 
of orbital d xy is adjusted from 0.3 to 0.38 and \x is adjusted 
to 0.006 in the upper panel. All energy quantities are in units 
of eV. 



appearance of the e-FS. To produce the e-FS, the on-site 
energy of orbital d xy is adjusted from 0.3 to 0.38, and 
H is adjusted from to 0.006 to keep the band filling at 
n = 6 in the upper panel, where the divergent suscepti- 
bility is X1414 at T = 0.019 and U = 1.38. In the lower 
panel, X1414 is suppressed by the annihilation of the e- 
FS, and the divergent channel changes to Xmi- With 
the definition of pseudo-transition-temperatures T% and 
T* 2 as X s (Q,T* d2) ) = 10 2 X°(Q,T* 1(d2) ), our numerical 



results show that the transition temperature T d2 = 0.019 
is higher than T dl = 0.012 in the presence of e-FS, which 
leads to a higher transition temperature of OTDW than 
that of SDW, thus the separation between Ts and Tjv 
appears, as is observed in 1111 compounds. Without e- 
FS, we get T dl = 0.017 > T* 2 = 0.0068 which indicates 
same T$ and Xjy in our scenario, as is observed in 122 
compounds. In this sense, the essential difference in the 
structure distortion between 1111 and 122 systems may 
be attributed to the existence of e-FS in 1111 parent 
compounds. Meanwhile, we want to emphasize that de- 
pending on the detailed band structure, there are other 
possible realizations of OTDW, such as the inter-orbital 
nesting between d xz and d yz orbitals. The importance of 
the e-FS in splitting Ts and Tjy awaits for further exper- 
imental verification. 



IV. SUMMARY 

We make several remarks as the summary. (1) As has 
already been addressed, the OTDW is not a static or- 
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dered state of total charge or spin density, which makes it 
"invisible" to the neutron detections. (2) The transverse 
nature of OTDW implies that the ordered orbital com- 
ponent is a mixture of the two involved orbitals, for in- 
stance, a mixture of d xz /d yz and d xy orbitals in the iron- 
based superconductors. (3) In our model calculations, 
the instabilities of OTDWs occur in the space spanned 
by d xz (or d yz ) and d xy orbitals. However, depending 
on the band structure, other possibilities exist. (4) The 
OTDW might be related with the "hidden" order claimed 
by some experiments (32|. (5) One intriguing feature 
stemming from the intinerant multi-orbital model is that 



the established SDW should be orbital-polarized (as is 
shown in Fig. 3). 
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